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We consider a model for a Planck scale ultraviolet cutoff which is based on Shannon sampling. 
Shannon sampling originated in information theory, where it expresses the equivalence of continuous 
and discrete representations of information. When applied to quantum field theory. Shannon sam¬ 
pling expresses a hard ultraviolet cutoff in the form of a bandlimitation. This introduces nonlocality 
at the cutoff scale in a way that is more subtle than a simple discretization of space: quantum fields 
can then be represented as either living on continuous space or, entirely equivalently, as living on any 
one lattice whose average spacing is sufficiently small. We explicitly calculate vacuum entanglement 
entropies in 1-1-1 dimension and we find a transition between logarithmic and linear scaling of the 
entropy, which is the expected H-l dimensional analog of the transition from an area to a volume 
law. We also use entanglement entropy and mutual information as measures to probe in detail the 
localizability of the field degrees of freedom. We find that, even though neither translation nor 
rotation invariance are broken, each field degree of freedom occupies an incompressible volume of 
space, indicating a finite information density. 

PACS numbers: 03.67.-a, 03.70.-l-k, 04.60.-m, 89.70.-a 


I. INTRODUCTION 

It is generally thought that quantum theory and gen¬ 
eral relativity, when combined, imply the existence of 
a minimum length in nature, see e.g., [TH6]. One line 
of argument is to consider that when attempting to re¬ 
solve distances with a smaller and smaller uncertainty 
Ax, the thereby necessarily increasing momentum un¬ 
certainty will induce a correspondingly increasing curva¬ 
ture uncertainty. Eventually the growing curvature un¬ 
certainty should make it impossible to resolve distances 
any further, leading to Ax > £p, where £p is at the 
Planck scale, or perhaps the string scale. 

Another line of argument comes from black hole 
physics. Black holes (and more general causal hori¬ 
zons H) are believed to carry a finite entropy given 
by the Bekenstein-Hawking formula, 5'bh = ^ (we set 
c = h = 1). This entropy may have its origin in the en¬ 
tanglement entropy of quantum fields [5HTU] (see [TT] for 
a review): the entanglement entropy has the same scaling 
behavior as ^bh and the dependence on the number of 
particle species could also be matched when taking into 
account how the renormalized gravitational coupling de¬ 
pends on the number of species [mils!. Crucially, the 
entanglement entropy also matches the order of magni- 


* j2pye@uwaterloo.ca 
^ donnelly@physics.ucsb.edu 
^ akempf@uwaterloo.ca 


tude of 5 bh; if there exists a natural ultraviolet cutoff 
which is close to the Planck length, £p = VG. 

Given these and other arguments, (see e.g., [H]), for 
the existence of a natural ultraviolet cutoff at the Planck 
length, the question arises as to the nature of this ultra¬ 
violet cutoff. In which sense might the density of degrees 
of freedom in quantum fields be finite? How are these 
degrees of freedom distributed and how can they be spa¬ 
tially resolved? 

To completely answer these questions would require 
knowledge of the still-to-be-developed theory of quantum 
gravity. More feasible at the present stage is to study 
how the natural ultraviolet cutoff may first manifest itself 
as one approaches the Planck scale from low energies, 
i.e., when coming from the safe ground of low energy 
physics where conventional quantum field theory (QFT) 
still holds. The question then is how can one model the 
first modifications to quantum field theory that express 
the impact of a Planck length cutoff? 

One simple model is that of QFT on discretized space 
or spacetime. A difficulty with this approach is the as¬ 
sociated complete breaking of local Poincare symmetry, 
though there are interesting methods to generate dis¬ 
cretizations whose statistics are Poincare invariant [H- 
HT]. Also, when lattices evolve in time, there tend to be 
problems with non-adiabaticity and an associated exces¬ 
sive particle production from the vacuum, in particular, 
in an expanding universe m- 

Here, we will consider a different model for a natu¬ 
ral ultraviolet cutoff in QFT with which at least the lo- 
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cal Euclidean symmetries are preserved and which does 
not necessarily induce non-adiabaticities in an expanding 
universe. The model is that of QFT with a hard cutoff 
in the form of a finite spatial bandwidth, i.e, a finite 
smallest wavelength How does this compare to 

the models of QFT on a lattice? It is well-known that 
the discretization of space gives rise to a minimum wave¬ 
length. It is less well-known that, vice versa, to impose a 
lower bound on wavelengths does not automatically lead 
to a discretization of space. Instead, one can also ob¬ 
tain a so-called sampling theoretic cutoff which does not 
break symmetries such as translation invariance. With 
this cutoff, the theory does not live on a single dedicated 
lattice. Instead, it lives on continuous space and, entirely 
equivalently, it can be written as living on any one lattice 
whose average spacing is sufficiently small. 

The mathematical structure that underlies the sam¬ 
pling theoretic ultraviolet cutoff is Shannon sampling 
theory. Shannon sampling plays a central role in in¬ 
formation theory because it establishes the equivalence 
of continuous and discrete representations of information 

Consider, for example, a music signal, /(t), that is low- 
pass filtered to some finite bandwidth, H. Even though 
the bandlimited music signal is a continuous function in 
time, it suffices to know the signal’s amplitudes f{tn) on 
any single lattice of points {tn}, if the lattice has an av¬ 
erage spacing at or below the critical spacing tt/H, which 
is called the Nyquist spacing. From the amplitude sam¬ 
ples {f{tn)}, the signal f{t) can then be reconstructed 
for all t, in principle, without error. This establishes the 
equivalence of continuous and discrete representations of 
information. The fact that none of the set of sufficiently 
densely-spaced lattices is preferred allows the preserva¬ 
tion of full translation invariance (and in higher dimen¬ 
sions also rotation invariance). 

Applied to physics, this means that spacetime could 
be simultaneously continuous and discrete in mathemat¬ 
ically the same way that information (such as a music 
signal) possesses simultaneously both a continuous repre¬ 
sentation and equivalent discrete representations. Phys¬ 
ical fields then possess a representation on continuous 
space, while being fully equivalently represented also by 
their amplitudes on any one lattice of sufficiently dense 
spacing. The democracy among these lattices allows 
translation and rotation invariance and more generally 
also Killing vector fields to be preserved with the cutoff. 

Our aim here is to study the implications of this sam¬ 
pling theoretic ultraviolet cutoff for localization and en¬ 
tanglement in quantum field theory. 

We find that the sampling-theoretic ultraviolet cutoff 
implies a particular kind of nonlocality in QFT which 
manifests itself through small but non-vanishing equal¬ 
time commutators at spacelike separation. This type 
of nonlocality appears naturally in perturbative quan¬ 
tum gravity [SUES] and in string field theory [26]. We 
also show how correlation functions are affected by the 
bandlimit, both classically at finite temperature, and 


quantum-mechanically. 

Further, in order to probe the spatial localization of de¬ 
grees of freedom with this UV cutoff, we examine the be¬ 
haviour of the entanglement entropy in a quantum field. 
For comparison, recall that in a QFT regulated on a lat¬ 
tice, the local field oscillators are coupled via the dis¬ 
cretized spatial Laplace operator in the Lagrangian. This 
makes their joint ground state, i.e., the vacuum state, an 
entangled state. Tracing over a region therefore yields 
an entanglement entropy. Due to the generally short 
range of the vacuum entanglement, the contributions to 
the entanglement entropy are dominated by correlations 
between those local field oscillators that are close the 
boundary of the considered region. This usually leads to 
an area law for the entanglement entropy, see e.g., m- 
In the 1-1-1 dimensional massless field theory that we will 
consider, the correlation length is infinite and the en¬ 
tanglement entropy therefore grows logarithmically (see 

m)- 

Unlike in such lattice theories, a bandlimited QFT does 
not have a preferred lattice representation. This makes 
the splitting of space into distinct regions nontrivial and, 
as a consequence, the Hilbert space does not automati¬ 
cally factorize into subsystems. Instead, more subtly, any 
discrete subset of points in space now defines a subsys¬ 
tem, and one can calculate its entanglement entropy with 
the rest of the system. In particular, we show how to cal¬ 
culate the entanglement entropy for a subset of degrees 
of freedom in a Gaussian state. Our result generalizes 
known formulae for the entropy of Gaussian states to the 
setting where the commutation relations are nontrivial. 

With samples placed at the Nyquist spacing, tt/H, we 
recover the usual logarithmic scaling behaviour of the en¬ 
tanglement entropy. When the spacing is larger than the 
Nyquist spacing (undersampling), we find that the entan¬ 
glement entropy crosses over to a volume law. Further, 
unlike in a lattice theory, in the bandlimited theory it is 
possible to probe the entanglement between degrees of 
freedom that are arbitrarily closely spaced. Surprisingly, 
when the spacing of samples is smaller than the Nyquist 
spacing (oversampling), we find no reduction in entan¬ 
glement entropy. Since the sampling points can occupy 
an arbitrarily small region, this, naively, appears to indi¬ 
cate that the field still carries an infinite density of local 
degrees of freedom (local field oscillators). We show that 
the resolution of this apparent paradox is that the region 
of space being probed by the sample points does not ac¬ 
tually decrease as the samples are taken closer together. 
In effect, each local field degree of freedom occupies an 
incompressible volume. 

Finally, we briefly examine the infrared behaviour of 
the entanglement entropy in the 1-1-1 dimensional ban¬ 
dlimited theory. In addition to the well-known logarith¬ 
mic growth of the entanglement entropy, there is a sub¬ 
leading double logarithmic infrared divergence. This di¬ 
vergence appears whenever there is a continuum of modes 
above the infrared cutoff scale; thus it appears when reg¬ 
ulating the infrared with a mass or a hard momentum 
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cutoff, but not if regulating by imposing periodic bound¬ 
ary conditions. 

We conclude with a discussion of implications and di¬ 
rections for future work. 


II. SAMPLING THEORY FOR QUANTUM 
FIELDS 

A. Overview of classical sampling theory 

The central result of classical sampling theory is Shan¬ 
non’s sampling theorem . This theorem establishes an 
equivalence between discrete and continuous representa¬ 
tions of information and it is, therefore, in ubiquitous 
use in communication engineering and signal processing. 
Let us consider a bandlimited signal, 4>{x), i.e., a sig¬ 
nal whose Fourier transform is supported on an interval 
(—Q, Q), whose width is the bandwidth 20: 

dk 

^ ( 1 ) 


not confined to any particular sampling lattice, the in¬ 
formation contained in bandlimited signal is subtly non¬ 
local. 

Sampling theory generalizes readily to bandlimited 
functions in higher dimensions, whose Fourier transforms 
are supported in a compact region of momentum space, 
see [SIH31 and it generalizes also to curved space [55] . 

Now we will briefly outline some of the functional ana¬ 
lytic properties of the space of bandlimited functions. We 
refer the reader to Appendixj^for more details regarding 
the functional analytic structure of sampling theory. 

The sampling theorem for a Nyquist lattice 
implies that the collection of sine functions centred at the 
lattice points forms a basis for the space of bandlimited 
functions. Moreover, the coefficients of a function in this 
basis are simply the values of the function at the corre¬ 
sponding lattice points. This basis is orthonormal in the 
inner product 

:=^ J dx (3) 

The orthogonality follows from the identity 


Shannon’s theorem states that such a function is com¬ 
pletely determined by its values on a discrete 

lattice of points {x^^}, where = (27rn — a)/(2fl) 
and n € Z. Here, a S [0, 27r) is an arbitrary hxed con¬ 
stant that labels lattices. Given the values, {(?!)(a:n“^)}, 
of the function on this lattice, (j){x) can be recovered for 
any x by the reconstruction formula: 


— J dx sinc[(x — sinc[(x — = Snm- (4) 

Different values of a corresponds to translated versions 
of the lattice. The sine functions centred at lattice points 
with a different value of a also form an orthogonal basis. 
Crucially, however, sine functions from different lattices 
are not orthogonal: 


<^(*) = H sinc[(x - (2) 


Thus, the space of bandlimited functions is completely 
equivalent to the space of functions defined on this lat¬ 
tice. This has the practical implication of allowing us 
to work with concrete and therefore computationally- 
convenient representations of functions on a lattice (as 
well as quantum fields, as we will see below) while pre¬ 
serving translation invariance (because the function re¬ 
tains its equivalent continuous representation). 

Another important finding in sampling theory is that 
the samples of the function do not have to be taken 
equidistantly: a bandlimited function can be recon¬ 
structed on any discrete set of points as long as these 
points are chosen with a sufficiently dense average spac¬ 
ing (technically, the Beurling average spacing [55115D]1. 
Here, this maximum average spacing is tt/Q. In gen¬ 
eral, such a set of points is called a sampling lattice, and 
the case of equidistant samples with separation tt/D is 
called a Nyquist lattice. In the case of sampling on a lat¬ 
tice other than a Nyquist lattice, perfect reconstruction is 
still possible but the reconstruction formula is more com¬ 
plicated than and the reconstruction becomes more 
sensitive to noise in the samples. Because the degrees of 
freedom of the function (values at the sample points) are 


n 

TT 



)Q] sinc[(a; 




= sinc[(4“^ 

= sinc[7r(n — m) + {a — Q!')/2] 
^ ^nm 5 fo^ OL ^ Oi 


(5) 


Here a, a' G [0,27r). This fact will be important when we 
study the localization of field degrees of freedom. 

Let us denote the vector in the function space corre¬ 
sponding to sinc[(a: — a;i“^)Q] by \x^'^). For hxed a, we 
then obtain a resolution of the identity: 

^|x(“))(:r(“)|=l (6) 


Taking the union of these bases over a, we get an over¬ 
complete basis for the function space with the corre¬ 
sponding ccn'^^’s covering R. Thus, we can write down 
an overcomplete resolution of identity: 


D 

TT 


dx |a;)(a;| = 1, 


(7) 


where the measure U.dx/'K gives the density of degrees 
of freedom in space, given by the Nyquist rate. This 
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overcomplete resolution of the identity is analogous to 
that for coherent states |a) 



where the measure Sol/'k = dxd’p/l'nh gives the density 
of independent states in phase space. 

The space of bandlimited functions is a reproducing 
kernel Hilbert space [3ll [32l |34j . This means that any 
function in the space can be recovered by means of the 
reproducing kernel K(x,x') as follows: 

= I 

K{x,x') := —{x\x') = — sinc[(a; — x')f2]. (9) 

TT TT 

B. Reconstruction formula for quantum fields 

We now apply the sampling theorem to a quantum 
field (j) in 1+1 dimensions, with a Hilbert space H. Let 
|(/)) S "H be an eigenstate of the field operator obeying 
(i){x) \(j)) = (j){x) \(j)) Va: e IR for some real-valued function 
of eigenvalues, tp. Note that in a continuum field theory 
without a cutoff, (p is not an operator but an operator¬ 
valued distribution which must be smeared with a suit¬ 
able test function to obtain an operator. In the ban¬ 
dlimited theory, the cutoff acts as an effective smearing 
function, and 0 is a genuine operator, though still un¬ 
bounded. 

Now let be the subspace spanned by the eigen¬ 

states of 0 where the corresponding functions of eigen¬ 
values 0 are functions bandlimited by H. Because 0 is 
a bandlimited function, the eigenvalue p{x) at any point 
a: € R can be determined by the knowledge of the eigen¬ 
values 0(a;^r^) at all of the points xlP^'^ on a sampling 
lattice {x^^}ne'z- Thus, the action of the operator 0(x) 
on an eigenstate of the field is determined from its action 
on a sampling lattice. Explicitly, 


p{x) 10) 

= (P{.X) 10) 




= E sinc[(a; 

- xy‘'>)n]p[xy‘y 10) 







= E sinc[(a; 

|0) . 

(10) 





Of course, this 

is true for all 

eigenstates of 0 in 'H( 

—o,n) ■ 


Since is the span of these eigenstates, the ac¬ 

tion of p{x) is determined by that of {0(a:n“^)}„gz for all 
states in Thus, we can write: 

sinc[(a: - a;£“))H]0(a;£“^). (11) 

Therefore, the operators {0(a:n“^)}„g^ form a complete 
set of commuting observables for any a. The fact that 


there are many lattices on which the field can be rep¬ 
resented means that the localization of the degrees of 
freedom is non-trivial. In the next section we will briefly 
examine this, but it will be studied further in later sec¬ 
tions of the paper. 


C. A first look at localization 

In a bandlimited field theory, the local harmonic os¬ 
cillators, i.e., the local degrees of freedom p{x), are not 
all independent. Namely, a set of degrees of freedom 
{p{xn)} is linearly independent only if the a;„ all belong 
to the same sampling lattice. The field p{x) at any other 
spatial point x can then be reconstructed from the am¬ 
plitudes {p{xn)}- In this section we will construct a a 
function that describes the spatial volume occupied by a 
set of degrees of freedom {p{x)}. 

To this end, consider the subspace spanned by N po¬ 
sition eigenvectors of first quantization (i.e., 

they can be represented as number-valued function over 
space). We will allow that they are not all from the same 
Nyquist lattice, and therefore they are generally not or¬ 
thogonal. Intuitively, the vector \xn) characterizes the 
spatial profile of the field degree of freedom 0(a;„) lo¬ 
cated at this point. Let us now construct the projector 
onto this subspace, in the basis of these N position vec¬ 
tors. First, we map the position vectors to an orthogonal 
basis 

\ei) Bij\xj). (12) 

j 

The projector onto this subspace is given by 
Iat := ^ |ei)(e^| 

i 

= E (13) 

j,k \ i / 

We can express the elements of the projector in the non- 
orthogonal basis {|a;„)}()0]^ using the reproducing kernel, 

O 

K{xj,Xk) ■■= —{xj\xk) 

TT 

= (14) 

TT 

i 

Viewing K(xj,Xk) as the {j,ky^ element of an iV x V 
matrix iLjv, we can write 

Iat = ~ E!(^Af )kj\xj){^k\- (15) 

TT ^' 

3,k 
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Inserting the resolution of identity Q, we can write 
this projector in the continuum basis: 
r ^ ^ 

In = dxdx' K{x,Xj){Kj})kjK{xk,x') |a;)(a;'| 

^ jk 

(16) 

Note that the trace of this operator is N, since it is sim¬ 
ply the identity on an TV-dimensional subspace. This 
trace is represented in the continuum basis as an integral 
over the diagonal elements of the integral kernel. We 
can interpret the diagonal elements of this projector in 
the continuum basis (i.e., as a function of a; S R) as the 
spatial profile of the N vectors {|a:„)}()Li- We can then 
visualise the spatial profile of N degrees of freedom for 
various spacings between the degrees of freedom; this is 
illustrated in Figure for 5 points. 


Ax 
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FIG. 1. Spatial profile (in the horizontal axis) for 5 degrees 
of freedom as a function of their spacing (vertical axis). Both 
axes are scaled so that the cutoff length, or Nyquist spacing, 
■k/Q, is equal to 1. The sample points are centred at a; = 
0. We see that for spacings above the Nyquist spacing, the 
degrees of freedom occupy a region consisting of 5 disjoint 
intervals of length ~ 1 surrounding each sample point. Below 
the Nyquist spacing, the degrees of freedom merge to occupy 
a single interval of length ~ 5. This interval does not decrease 
in size even as the sampling points are taken on top of one 
another. 


We see that for sampling point separations larger than 
the Nyquist spacing (shown in the region > 1 on the ver¬ 
tical axis), the degrees of freedom are clearly localized 
about the chosen sampling points and they are nonlocal 
on the order of the cutoff length (in the horizontal axis). 
When the points are spaced close to the Nyquist spac¬ 
ing (i.e., spacing of 1), the spatial profiles of the points 
converge to form a volume of length 5. The surprising 
point illustrated by this plot is that when the points are 
spaced closer than a Nyquist spacing apart, the points 
still describe a spatial volume equivalent to the spatial 
volume described by the points when they are spaced at 
the Nyquist spacing. Thus we see that each point can 


be localized to a unit of volume, and that these volume 
elements are incompressible, in the sense that N points 
pushed close together still describe an interval of length 
TV. Intuitively, it is clear that, because of the cutoff, 
placing the sampling points closer than a Nyquist spac¬ 
ing should not uncover new degrees of freedom. What is 
surprising but in hindsight plausible is that even when TV 
sample points are closer than the Nyquist spacing they 
still access TV degrees of freedom. 

We will return to this issue in the context of quantum 
fields in section IVl 


D. Bandlimited correlation functions 


We would now like to see the effect that bandlimiting a 
quantum field has on the commutation relations and cor¬ 
relation functions of the ground state of the field. The 
particular fields we will hereafter be considering are 1-1-1 
dimensional massless scalar quantum fields with an ultra¬ 
violet cutoff n and infrared cutoff uj imposed on the spa¬ 
tial momentum, so that w < |fc| < R. Note that bandlim¬ 
ited functions which have support on (—fl, —cu) U (uj, H) 
in Fourier space form a subspace of the space of func¬ 
tions with support on (—fl, 11). Thus the above sampling 
theorem § also applies to these functions. 

We consider the system described by the free Klein- 
Gordon Hamiltonian, 


f dx (■n'^(x) + (t)(x)(-A)(t)(x)) , (17) 

^ Jtr. 


where A = 9^ is the scalar Laplacian operator. Here and 
below we will not write the " on top of operators. The 
field (j) and its conjugate momentum tt can be expressed 
with the usual mode expansion (see, e.g., [3H]); except 
that the field modes outside of the range oj < |fc| <11 are 
removed: 


Hx) 

Tr(x) 


dk 1 / 


L<\k\<Q. 27r ^/2\k\ 
f dk 1 
luj<\k\<fi 27r i 


[ake 


ikx 


( 


Q'kC 


ikx 


+ , (18) 

-4e-*^“).(19) 


Here ak,a\ usual annihilation and creation op¬ 

erators obeying the commutation relations [ak,a\,] = 
(27r)S(k — k') and all other commutators vanish. No¬ 
tice that the Hilbert space is unitarily preserved 

under free time evolution because the Fourier modes of 
the field are uncouplec^^ From the mode expansions we 


^ In an interacting theory the modes would couple and, naively, one 
may expect this to generate shorter-than-cutoff (i.e. transplanck- 
ian) wavelengths. Realistically, however, such high-energetic par¬ 
ticle collisions would necessarily excite the very gravitational de¬ 
grees of freedom that are thought to enforce the ultraviolet cutoff 
in the first place, which may well save unitarity. 
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can now calculate commutation relations and two point 
functions for the field at two arbitrary points x, x' € R. 
The equal-time commutation relations between (f> and tt 
become 

), ( 20 ) 


[4>{x),tt{x')] = i ( — sinc(nAa;) — — sinc(a;Aa;) 

\Tr TT 


where we have written Ax := x — x'. Note that the com¬ 
mutator can be related to the reproducing kernel which 
we used in the first quantized picture ([^ after projecting 
the kernel onto the space of bandlimited functions with 
support on (— —uj) U (w, fl) in Fourier space: 

[(j){x)^TT{x')\ = iK{x,x). (21) 


The remaining commutators all vanish: [(j){x),(j){x')] = 
[7r(x), 7r(a:')] = 0 for any two points. 

The correlation function for the ground state of the 
field at two distinct points is: 


(0| (j){x)(j){x') |0) = — [Ci(r2Ax) — Ci(a;Aa:)]. 

ZTT 


In the coincidence limit this becomes: 


l<(’^{a;) |0) = ^log 


( 22 ) 


(23) 


Similarly, the correlation function for the conjugate mo¬ 
mentum at two distinct points is 


17r(a:)7r(a:') |0) = 


cos(ilAa:) — cos(a;Aa:) 
27rAx^ 

n sin(r2Aa;) — w sin(a;Aa;) 


(24) 


2ttAx 

which becomes in the coincidence limit 


K2(x)|o) = 




47r 


(25) 


The correlation functions between cj) and tt at equal times 
all vanish, ^ (0| {^(a:), 7r(x')} |0) = 0. 

From Eq. (231, we see that the infrared cutoff w is nec¬ 
essary to regulate the infrared divergence of the cor¬ 
relator. It will therefore be necessary to keep uj strictly 
positive for our calculations, and we will examine infrared 
effects in more detail in section IVII Since we are more 
interested in the effect of the ultraviolet cutoff fl, we dis¬ 
play the correlation function in the limit w —>■ 0 while 
adding appropriate counterterms to the correlation func¬ 
tions. The second term of Eq. (22) contains the infrared 
divergence. For w/fl I, we have 

Ci ^^nAx^ = 7 -|-log(nAa;)-|-log(a;/n)-|- 0 (a;/n). (26) 

Therefore if we add a term log(w/f2)/(27r) to the {(fxj)) 
correlator, in the limit a; —>■ 0 we find that at two distinct 
points 

(0| (j>{x)(j){x') |0) I—)■ — [Ci(nAa;) — 7 — log(nAa;)]. (27) 
27r 


We also want this correlator to give the usual nonban- 
dlimited correlator when fl —> 00 . The nonbandlimited 
correlator is 


r rjh ^ . _i 

(0| ^(a:)(/)(a;') |0) = ^ + log |Ax|). 

(28) 

Thus, since Ci(n|Aa:|) —>■ 0 as 11 —>■ 00 , we must add 
another factor of log(n)/(27r) to (</)(/>) to get 

(0| (j){x)(l){x') |0) = ^[Ci(IIAa:) - 7 - log |Aa;|], (29) 

and at a single point we have 

(0|<()2(a:)|0) = ^log(II). (30) 

The correlator {tttt) is not divergent in the infrared, thus 
we can directly take the limit w/fl —)■ 0 . 

In Figures and we present graphs of the correla¬ 
tion functions (0| (()(a;)(/i(a;') |0) and (0| 7 r(x) 7 r(a;') |0) (re¬ 
spectively) as a function of Ax, with the infrared cutoff 
removed. In both figures, we see that the correlations 
decay with distance, similarly to the correlations when 
there is no ultraviolet cutoff. We notice that the ban¬ 
dlimited correlation functions oscillate with wavelengths 
of the order of the ultraviolet cutoff length. As we will 
see later, similar oscillations also occur in the correlations 
of thermally distributed bandlimited classical signals. 



FIG. 2. correlations as a function of their separation. 

The horizontal axis is scaled by fl/rr so that integer values 
correspond to Nyquist spacings. The bandlimited correlations 
are blue with the Nyquist spacings indicated by red dots. The 
black dashed line shows the ultraviolet-divergent correlations 
without the ultraviolet bandlimit. We see that for points on 
the Nyquist lattice, the bandlimited correlators are in closer 
agreement to the correlation functions without the bandlimit. 
A counterterm of log(a;) / (27r) is added to {(jxj)) to cancel the 
infrared divergence. 
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FIG. 3. TT-TT correlations as a function of their separation. 
The horizontal axis is scaled by SI/tt so that integer values 
correspond to Nyquist spacings. The vertical axis is scaled by 
The bandlimited correlations are blue with the Nyquist 
spacings indicated by red dots. The black dashed line shows 
the ultraviolet-divergent correlations without the ultraviolet 
bandlimit. 


III. ENTROPY OF GAUSSIAN STATES 


We would like to understand how the delocalization 
of degrees of freedom of the bandlimited quantum field 
theory affects the localization of information in space. A 
useful tool to probe the distribution of degrees of freedom 
in space is the von Neumann entropy of localized subsys¬ 
tems, and derived quantities such as the mutual infor¬ 
mation. Suppose we have a state of a bandlimited field, 
whose wavenumbers are in the range a; < |fc| < fl, and 
we are able to make measurements of the field amplitude 
(j) and its conjugate momentum tt at a finite number, TV, 
of points. Any set of points defines a subsystem in this 
way, regardless of their relative positions in space, and 
we can define a von Neumann entropy associated with 
this subsystem. In this section we show how to calculate 
the von Neumann entropy associated with an arbitrary 
subset of points using the formalism of Gaussian states 
1571 . We leave concrete numerical calculations for section 

m 


First in section [III A[ we consider a classical field in a 
thermal state. For this system there is no entanglement 
or quantum noise and all of the entropy comes from ther¬ 
mal fluctuations. Then in III B we consider a thermal 
state of a quantum system, which in the limit of zero 
temperature becomes the ground state of the quantum 
field. In this case the entropy includes both entangle¬ 
ment entropy and thermal fluctuations, and in the zero 
temperature limit is entirely due to quantum entangle¬ 
ment. This will demonstrate how the classical results 
emerge from the quantum results at high temperatures. 


A. Classical entropy 


Let us first consider the classical situation where the 
field we are given is chosen from some given classical 
distribution. If, for example, the statistical distribution 
is that of Gaussian white noise, the values of the field 
at different sample points are uncorrelated. Since the 
entropy is additive for uncorrelated degrees of freedom, 
the entropy associated with N samples is linear in N. 
The classical case which we will be more interested in, 
because it is more comparable to the quantum case, is 
the case where the statistical distribution of the signals 
is a thermal distribution of a 1-1-1 dimensional classical 
Klein-Gordon Hamiltonian, with fields that have both 
an ultraviolet cutoff H and infrared cutoff w. For this 
distribution, the degrees of freedom of the field at the 
sample points will be coupled, causing the values of the 
field at these points to be correlated. We therefore ex¬ 
pect that the entropy of the interval of N samples has a 
contribution which is nonlinearly dependent on N. We 
now explicitly calculate the entropy created for these N 
samples. 

We will go through the entropy calculation in detail 
since many of the steps are reproduced in the quantum 
calculation. We begin with the Hamiltonian in momen¬ 
tum space, 

= l f ^ + • ( 31 ) 

^ Jij<\k\<n V / 


Both the ultraviolet and infrared cutoffs are explicitly 
enforced in this description by the restriction on the mo¬ 
mentum k. For a thermal state, the probability distri¬ 
bution of the fields in momentum space is given by the 
Boltzmann distribution 

= (32) 


We seek the entropy of a reduced probability distribu¬ 
tion for a set of N sample points. To this end, we choose 
a lattice {a;„}„g^ and we write Axmn ■= Xm — Xn- We 
denote the values of the field at these points (j)n ■= 

TTn ■= 7r(xn)- To find the probability distribution for the 
samples, we perform a change of phase space coordinates 
from the fields in momentum space to a sampling lattice: 


^(Xn) 

^(Xn) 


> uj<\k\ <n 

/ 

' ci;<|fc| <n 


fih 


(33) 


Instead of explicitly performing this change of variables, 
it will be more convenient to make use of the fact that the 
probability distribution in terms of (j)n and 7r„ is Gaus¬ 
sian. This follows from the fact that the probability dis¬ 
tribution in momentum space is Gaussian and the change 
of phase space variables (33) is linear. The distribution 
is therefore characterized entirely by the two-point func¬ 
tions. 











First, note that one can easily calculate the power spec¬ 
tra of the fields from the probability distribution (32) 
where the fields are represented in momentum space: 


mm = (34) 

p^k 

mm = (35) 


Then we can obtain the correlators of the fields between 
two arbitrary points by taking the Fourier transform of 
the power spectra: 




cos(a;Aa;) cos(f^Aa:) 

TTU} TtVI 


/\ nr 

[Si(a;Aa;) — Si(nAa;)], 


(36) 


{'n{x)Tl{x')) = - 


^ A \ ^ / A \ 

— sinc(i 2 Aa:)-smc(cjAa:j 

TT TT 


•(37) 


These correlators are plotted as a function of Aa; in 
Figures and alongside the finite-temperature quan¬ 
tum correlators calculated below. In a similar manner 
as for the vacuum correlators in the previous section, 
these functions are displayed with added infrared coun¬ 
terterms. The (tttt) correlator is again finite as w/n —>■ 0, 
so for plotting purposes, we can just take this limit. Since 
Si(a;Aa;) —>• 0 as w/fl —>■ 0, the first term in Eq. ( [M| ) is 
the only infrared divergent term, whose divergence may 
be cancelled by subtracting since as w/fl —>■ 0 , 


cos(a;Aa:) 


TT/Sw TTj3uj 

The function plotted in figure is then 


imm)) = 


cos(nAa:) 

Tril 


Aa; Si(nAa;) 


(38) 


(39) 


For the entropy calculation, we require the reduced 
probability distribution for N sample points obtained by 
marginalizing over a complementary set of phase space 
coordinates. In a lattice theory, this would simply mean 
integrating over the field values at all lattice sites not 
under consideration. In a bandlimited theory, the com¬ 
plementary subsystem to a set of points is not generally 
associated to any set of points; it is nonlocal. As pre¬ 
viously noted, we do not actually need to carry out the 
marginalization, as the reduced probability distribution 
on the N samples is Gaussian, and hence entirely deter¬ 
mined by the two-point functions of the N points. Nev¬ 
ertheless it is instructive to carry out the decomposition 
of the total system into the N sample points and their 
complement. 

We will perform this calculation for any Gaussian state 
of the form 

Pimm) = 

Z (40) 
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with positive-definite, symmetric matrices A and B. For 
the above example, the ((in’s and 7r„’s would represent 
the field amplitudes on a sampling lattice. We can split 
the phase space into a subsystem describing the N sample 
points and its complementary subsystem by performing a 
change of variables which splits the matrices encoding the 
Poisson bracket and symplectic form into a direct sum of 
the matrices which act on the individual subspaces sep¬ 
arately. It will be convenient to write the phase space 
variables in a vector r = (((i, tti, (( 2 , 7 r 2 ,...). Then the 
Poisson bracket for the total phase space with the ul¬ 
traviolet and infrared cutoff constraints enforced can be 
encoded in an anti-symmetric matrix Kij := {vi^rj^pB^ 
where 

{0i,7r,}pB = — sinc(nAa;i,)-sinc(a;Aa;io), 

TT TT 

{((ij, 4’j}pB = {TTi, tTjIpb = 0. (41) 


Note that this consistent with the quantum mechanical 
commutation relations ( 20 ). 

First we write the Poisson bracket as a block matrix. 


A = 


a T] 
T 

-77 7 


(42) 


where = —a and 7 ^ = — 7 , ensuring that A^ = —A. 
The blocks are arranged such that the upper left block 
contains the first 2N indices which correspond to the 
JV sample points which will remain after marginalizing. 
Now consider a change of phase space variables r' = Qf, 
where 


Q:= 



0 

I 


(43) 


The transformed Poisson bracket is 


A' := QAQ^ 


a 0 
0 r]’^a~^r] + j 


(44) 


Now the phase space splits as a direct sum of the N 
degrees of freedom remaining after marginalizing and its 
complement, with symplectic form 


N 

F .— ^ ) ((U )nmd(/i^ A 

n,m—l 

+ [(??'^o"N + 7 )"^]nmd</>'„ Ad7r;. 

(45) 


This construction requires that the matrix a be invert¬ 
ible, but we will now show that in the bandlimited theory 
this is always the case. We can express a as the matrix 
Oy = K{xi,Xj) = ^{xi\xj), where \xi), \xj) are the ban¬ 
dlimited position eigenfunctions defined in |II C| Suppose 
that a is not invertible. Then we can find a vector ai 
such that 

0 = aia^jQj = '^{x^\aiaj\xj) = || ^ai|a;i)|p, (46) 
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and hence ai\xi) = 0. This implies that for any ban- 
dlimited function /, we have 

0 = (/I ai\xi) = Y = Y (47) 

i i i 


However given any finite set of points Xi and target val¬ 
ues Ui, we can find a bandlimited function / such that 
f{xi) = Ui [35] ■ This contradicts (47), hence a must be 


invertible. 

Now, the probability distribution after performing the 
marginalization is of the form: 


7*red((^^’ i) 


^red 


X e 2 




(48) 


where H|jv and B\ff are positive-definite symmetric ma¬ 
trices given by 


(■BI )mn = (^m </>«), (49) 

(^IaT )'mn — • (^^) 


These correlators are just ( |36| ) and ( |37| ) evaluated at the 
points X = Xjn and x' = a;„. Note that the matrices 


and BjjY are simply the original matrices A and B with 
indices restricted to the subsystem m,n G {!,... ,7V}. 

Before we can calculate the partition function by in¬ 
tegrating over the phase space variables {4>n,'^n}n=i, we 
must first choose a measure for this phase space. When 
calculating the entropy for a classical probability distri¬ 
bution with continuous random variables, there is an am¬ 
biguity in the entropy caused by an ambiguity in the 
choice of measure. Here, we will fix the measure to 
be proportional to the symplectic form in this subspace 
(Eq. ( |4^ ) d(j)n A d7r„/(27r). This is sufficient to resolve 
the ambiguity in the entropy, and ensures that the en¬ 
tropy invariant under symplectic transformations of the 
phase space. As we shall see, the factor of 27r gives an 
entropy which matches the high temperature limit of the 
von Neumann entropy in the quantum setting. 

Now we find that the reduced partition function is 


ydet(a) [[ ^ 


= det(a) ^\/det (A“i) det 


(51) 


Then the entropy of the distribution of the remaining 
samples is: 


N 


S = I det(a) ^ Jl 7^red({<(’"’’^"}"=i) 


n—1 


X log [Pred(('(’"’^"J-n=i)] 

= iV-I-log det(a)“^-\/det (A“i) det {B~^) . (52) 


Note that the first term is linear in the number of sam¬ 
ples TV, but the second term, which encodes the coupling 


between the sample points, is not exactly linear. In Fig¬ 
ure [^ we plot the entropy for a massless bandlimited 
field in a thermal state. We find that the entropy is very 
close to linear in the number of lattice points, regardless 
of the temperature. This simply reflects the extensivity 
of thermal entropy. 



FIG. 4. Entropy of a set of sample points for a thermally- 
distributed classical bandlimited Klein-Gordon field for sev¬ 
eral temperatures. We see that the magnitude of the entropy 
grows linearly with the number of points, illustrating a vol¬ 
ume law. Notice also that this figure shows negative entropy 
at low temperatures, which simply reflects the fact that con¬ 
tinuous probability distributions can have negative entropy. 
We fix the inherent ambiguity in these entropies by requiring 
that it matches the quantum entropy at high temperature. 


It turns out to be more convenient to rewrite the for¬ 
mula for the entropy in terms of the following matrices: 


AIat := 

SU := 


0 A 
-A 0 


0 A|j 


BIVA 0 
\n 


(53) 

(54) 


where [at denotes restriction of the matrix indices to 
1,... ,TV (i.e. the indices corresponding to the remain¬ 
ing sample points after marginalizing). The matrix 
encodes the Poisson bracket with the indices of the phase 
space vector Prearranged so that all of the (j)nS occur be¬ 
fore the TTji’s, and where Xmn '■= T^n}pB- Note that 
A is just the matrix a from Eq. (42) after a permutation 


of the indices. Similarly, the matrix T,\n is simply the 
matrix of correlators between the phase space variables. 

It is then straightforward to verify that the entropy 
can be written as 


S = N + log 

N 


\/det(A\YxY) 

= ^(I + log(d,)). 


(55) 
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where {di}^i are the positive imaginary parts of the 
eigenvalues of A|^^E|jv which come in pairs ±idi. There¬ 
fore, the entropy is determined entirely in terms of the 
Poisson bracket (41) and two-point correlators (36) and 
(37). Notice that since we can write down the bracket 
and correlators between any two points of the field (not 
necessarily on a Nyquist lattice), we can calculate the 
entropy associated with N arbitrary points of the field. 


B. Quantum entropy 


Let us now proceed to the case of calculating the corre¬ 
sponding entropy for N samples of a quantum field. Here 
we will review a method for calculating the entanglement 
entropy for Gaussian states of systems of quantum har¬ 
monic oscillators. The Gaussian formalism for calculat¬ 
ing the entropy for Gaussian states goes back to iHin] 
with significant simplifications developed in [371l39j . Gru- 
cial for our purposes is that this formalism also general¬ 
izes naturally to the case where the commutation rela¬ 
tions between (j) and n are nonlocal. As in the classical 
case, the state of a region, and hence its entropy, is de¬ 
termined entirely by the matrix of Poisson brackets and 
the covariance matrix. 

We begin with the Hamiltonian for a set of M coupled 
oscillators (where M may be finite or infinite): 



M 

E 




m-^mn n 


M 

E 




mn'vn • 


(56) 


Here A and B are positive-definite, symmetric matrices. 
As in section |HI A| it will be convenient to combine (f>n 
and TTn into a single vector f= {(pi, tti, (f> 2 ,'^ 2 , ■ ■ ■)• Then 
we define the correlation matrix 


■■= ^{{ri,rj}) = Itr {{ri,rj}p), (57) 

where p is the density matrix for the state. Also, since 
the commutators are c-numbers, we can define the matrix 
A that encodes the commutation relations: 




(58) 


simply factor into a tensor product of Hilbert spaces gen¬ 
erated by the field operators at the sample points; i.e., 
^(-n,n) 7^ This is because for points, say Xn and 

Xm, which do not lie on the same Nyquist lattice, the op¬ 
erator (p{xn) acts non-trivially on since this sector 
of the Hilbert space is generated by <p{xm) and 7r(xm), 
but [(j){xn),'x{xm)] 7 ^ 0. However, if we identify a finite 
subset of lattice points it is possible to factor the Hilbert 
space into a tensor product of a subspace describing the 
subsystem and its complement. 

Factoring the Hilbert space is similar to the phase 
space splitting performed in the classical case. We pro¬ 
ceed using the change of variables defined by the matrix 
(43), so that A takes the form: 


A' := QAQ^ = 


a 0 

0 r]’^a~^r] + j 


(59) 


Note that the block a := A|j corresponding to the 
degrees of freedom is unchanged under this 

transformation. The covariance matrix also transforms 
as 


E' := QEQ^. (60) 

It is also easy to check that the block E|j is unchanged 
under this transformation. 

Now, by Darboux’s theorem, we can find a change of 
variables via a transformation T which brings the matrix 
A' into its canonical form. That is, we can find a set of 
canonical coordinates satisfying canonical commutation 
relations [qi,pj\ = iSij. Since we have split the phase 
space into two pieces, we can find such a matrix T in the 
form 


Now, 


Ti 0 
0 Ta 


A" := TA'T^ = 


E" := TE'T^ 


0 


0 1 
-1 0 


(61) 

(62) 

(63) 


Now we would like to find the state of a subset of 
oscillators, i = 1,... ,N (with N < M) obtained by 
tracing over the complementary set of degrees of free¬ 
dom. For the bandlimited theory, the complementary 
set of degrees of freedom are not simply the complemen¬ 
tary set of lattce points N -I- 1,... ,M. Notice that if 
we remove the infrared cutoff, i.e., set w = 0, then on 
a Nyquist lattice we have the canonical commutation 
relations [(l){xm),'!^{xn)] = Therefore, on this 

Nyquist lattice the Hilbert space factors into the 

tensor product ®nB.x„ of the Hilbert spaces gener¬ 
ated by the field operators {^(a:„),7r(a:„)}. However, for 
a general sampling lattice the commutation relations are 
nonlocal at the cutoff scale. Thus, in the case of non- 
Nyquist sampling, the Hilbert space does not 


An advantage of working with covariance matrices is that 
the tracing operation can be implemented by simply re¬ 
stricting the indices of A" and E" to the set {!,..., 2A^}. 
Also, by construction of the matrices Q and T, the re¬ 
duced commutator and covariance matrices can be deter¬ 
mined solely from the reduced commutator and covari¬ 
ance matrices of the original variables, i.e., 

A'V = TiA|^Tf (64) 

E"U = TiE|ArTf. (65) 

Therefore, we see that since we are only mixing the coor¬ 
dinates for the degrees of freedom 1,..., A^ among them¬ 
selves before the tracing operation, the Hilbert space af¬ 
ter the tracing operation is the same as the Hilbert space 
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generated by the original operators {4'n,'^n}n=i- The 
importance of this fact is that it does not matter where 
the remainder of the samples are taken, as we can simply 
identify the entire Hilbert space Ets the tensor 

product of the Hilbert space generated by the operators 
corresponding to the Hilbert space of 

the complementary set of degrees of freedom. This also 
allows us to choose a set of samples which do not all lie on 
a Nyquist lattice, as we can simply identify the reduced 
Hilbert space as the space generated by the operators at 
these points. 

Now that we have the reduced commutator and co- 
variance matrices in a Darboux basis, as shown by 
Williamson [40l [41], we can make a further (symplec- 
tic) transformation S that preserves the form of A"|jv, 
and puts S"|Ar into diagonal form: 


SA"\nS^ 

SY."\nS'^ 


A" 


w, 


0 di 0 

0 d. 


( 66 ) 

(67) 


Here the diagonal entries di come in pairs, and are called 
the symplectic eigenvalues of S"|Ar. 

In these coordinates, the density matrix is a prod¬ 
uct of uncorrelated thermal density matrices of N har¬ 
monic oscillators with canonical coordinates pi satis¬ 
fying canonical commutation relations, in a state where 
(ft) = {Pi) = di- In particular, the oscillator is in 
the state 


= 


n>0 



|n)(n|. 


( 68 ) 


Note that the uncertainty relation implies 

di = AqiApi > l\[qi,Pi]\ = ^. (69) 


Thus the symplectic eigenvalues are all bounded below 
by di > i. 

The entropy is then the sum of entropies of each indi¬ 
vidual oscillator: 


s = Y,s{d,), 

i 

S{d) := {d + D log(d + D - {d 


(70) 

i)log(d- i). 


Thus we can determine the entropy entirely from the 
symplectic eigenvalues di. 

Note that one does not need to carry out this symplec¬ 
tic diagonalization in order to find the symplectic eigen¬ 
values. Under the sequence of similarity transformations 
Q, T, S, the eigenvalues of A|)^^S|jv are invariant, i.e.. 


spec 



di 0 
0 di 


= spec {[iST,)A\N{STif]-\ST,)i:\NiST,f) 

= spec (AD^^EIat) , 

( 71 ) 


where spec(A) denotes the spectrum of A. We therefore 
see that the eigenvalues of A|)^^S|Ar coincide with the 
eigenvalues of the matrix 


A'U' 


N 

0 

i=l 


di 0 
0 di 


0 


i=l 


0 di 

-d, 0 ’ 


(72) 


which are the values Aidi. Thus we can find the sym¬ 
plectic eigenvalues simply by hnding the eigenvalues of 
A|)^^E| 7 v- See [37], as well as the related result [33] . 

When the symplectic eigenvalues {di} are large (which 
in a thermal state can be shown to correspond to a 
high temperature limit for thermal states), the expres¬ 
sion S{d) becomes 

5(d)«l + log(d), (73) 


so the quantum entropy formula approaches the classical 
formula ( jS^ obtained above. This fact is illustrated in 
figure]^ which shows the classical and quantum formulas 
for S{d). 



d 


FIG. 5. The entropy of a harmonic oscillator as a function 
of the symplectic eigenvalue, both classically and quantum- 
mechanically. In quantum mechanics, the uncertainty prin¬ 
ciple requires d > 1/2, which is saturated by the vacuum 
state for which the entropy is zero. Classically, the symplec¬ 
tic eigenvalue can be any positive number, but the entropy 
becomes negative for small d. This again is a consequence of 
the fact that the entropy of a continuous probability distri¬ 
bution is not bounded from below. 


IV. VON NEUMANN ENTROPY 

We now calculate the entropy associated to a uniformly 
spaced lattice of field samples. There are three distinct 
regimes, depending on whether the spacing is equal to the 
Nyquist spacing, larger than Nyquist (undersampling), or 
smaller than Nyquist (oversampling). We will consider 
each of these possibilities in turn, both in vacuum and at 
finite temperature. 
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A. Nyquist sampling 

As shown in section the entropy of a set of sample 
points in a Gaussian state is determined by the commu¬ 
tators and two-point functions at the sample points. We 
now calculate these for a quantum held at hnite temper¬ 
ature. 

Before we perform the tracing operation, the density 
matrix associated with a single mode of the held is 


^ LXJ 


nfc=0 


Formally, the total density matrix for the held is 
P= 0 Pk- 

This leads to the power spectra 


(|<^feP)=tr(|</.fcpp) = - 
(kfeP) = tr(|7rfcpp) = Wfe 


1 


1 

-K - 


Wfe - 1 2y ’ 

1 r 

6/5“*= - 1 2, 


(74) 

(75) 

(76) 

(77) 


From the power spectra we hnd the two-point functions 
for a massless bandlimited held are 


^ +1] (78) 


=/3fe _ 1 2 


fji. 1 

{(j){x)(t){x'))= — cos(fcAa;)- 

Ju> k 

(7r(x)7r(x')) = J * cos{kAx)k ^ -(79) 

These correlators are plotted in Figures and [^alonside 
the classical correlators. We notice that, as expected, 
the classical and quantum correlators agree at high tem¬ 
peratures. Both classical and quantum correlations ex¬ 
hibit oscillations at the ultraviolet scale. At temperatures 
T/kl < 1, the quantum held exhibits stronger correlations 
than its classical counterpart due to vacuum huctuations. 

Using these thermal correlation functions together 
with the commutation relations (eq. (20)), we numeri¬ 


cally calculate the entanglement entropy of a hnite lattice 
of points separated by the Nyquist spacing tt/U. Then 
we consider the change in the entropy as the number of 
points increases. The entropy as a function of the num¬ 
ber of sample points is illustrated in Figure for a range 
of temperatures. The plot shows clearly the transition 
between linear and logarithmic growth of the entropy. 

Comparing the entropy in the quantum theory (Fig.j^ 
with the corresponding classical result (Fig. |^, we see 
the expected agreement of the scaling behaviour at high 
temperatures, with the entropy growing linearly with the 
number of lattice points. This simply reflects the exten- 
sivity of the thermal entropy. However, at low tempera¬ 
tures the entropy of the quantum field increases only log¬ 
arithmically with the number of points removed (e.g. the 
zero temperature state in Figure]^. This is a well-known 



FIG. 6. Classical and quantum correlations as a function 
of their separation at various temperatures. The horizontal 
axis is scaled by U/tt so that integer values correspond to 
Nyquist spacings. The vertical axis is scaled by U/T. The 
classical correlator is plotted as black dots. The scaling of 
the vertical axis absorbs all of the temperature dependence 
of the classical correlator (since it grows proportionally to 
T/U), thus the single graph of the classical correlator com¬ 
pletely characterises its behaviour. The quantum correlators 
are shown as lines for temperatures up to T/U = 1, at which 
point the quantum correlator converges to the classical corre¬ 
lator. 


result for conformal field theories, such as the massless 
scalar field, for which the entropy scales as S' ~ ^ log(UL) 
|28) . Here, for the thermal state at zero temperature (i.e., 
the ground state) the htted curve to the data is 

S{N) = 0.334 log A -h 3.25 (80) 

which is in agreement with the expected leading order 
term |log(AI). This result is for an infrared scale of 
w/H = 10-3OO. 

We therefore see that the total entropy of the quantum 
field is a combination of the thermal and entanglement 
entropy. At high temperatures the entropy is primarily 
thermal, while at low temperatures the thermal entropy 
vanishes, leaving only the entanglement entropy. 


B. Undersampling 

Now we examine how the entanglement entropy of the 
ground state of the quantum field scales with the number 
of samples N when the distance between adjacent sam¬ 
ples is independent of H. In sampling theory terminology, 
this corresponds to undersampling when Ax > tt/H and 
oversampling when Ax < tt/H. 

First, we will examine the case of undersampling. Re¬ 
call that the strength of correlations between the field at 
two separate points decays with the distance between the 
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FIG. 7. Classical and quantum tt-tt correlations as a function 
of their separation at various temperatures. The horizontal 
axis is scaled by SI/tt so that integer values correspond to 
Nyquist spacings. The vertical axis is scaled by Il/T. The 
classical correlator is plotted as black dots. The scaling of 
the vertical axis absorbs all of the temperature dependence 
of the classical correlator (since it grows proportionally to 
T/n), thus the single graph of the classical correlator com¬ 
pletely characterises its behaviour. The quantum correlators 
are shown as lines for temperatures approaching T/fl = 1, at 
which point the quantum correlator converges to the classical 
correlator. 


points and the entropy is dominated by local correlations. 
Thus, for points separated by more than the Nyquist 
spacing any given point is most strongly correlated with 
degrees of freedom in the complementary subsystem. As 
a result, the entanglement entropy in this regime is pro¬ 
portional to the number of points traced out. Hence we 
should recover a volume law for the entropy similar to 
the high temperature thermal state. 

The procedure to calculate the entanglement entropy is 
the same as before, except now the correlators and com¬ 
mutators are taken between points which do not all lie on 
a Nyquist lattice. Figure]^ shows the dependence of the 
entanglement entropy on the number of points traced out 
for spacings that vary between 1 and 2 Nyquist spacings. 

We see that for Nyquist and near-Nyquist sample spac¬ 
ings the entanglement entropy grows logarithmically with 
the number of points. As the sample separation increases 
there is a transition to a linear scaling with the number 
of points. This corroborates our intuition for the scaling 
behaviour of the entanglement entropy when undersam- 
pling. 

C. Oversampling 

Instead of separating the sample points, we now take N 
samples which are spaced closer than a Nyquist spacing 
apart. Starting from Nyquist spacing and pushing the 



FIG. 8. Entropy for N Nyquist spaced sample points for 
a thermal state of a quantum Klein-Gordon field. The plot 
illustrates the transition between the logarithmic behaviour at 
low temperature to linear behaviour at higher temperatures 
where the entropy becomes an extensive variable. 



FIG. 9. The dependence of the entanglement entropy on the 
number of sample points, for sample spacings between 1 and 
2 Nyquist spacings. The Ax labels in the legend are scaled by 
H/tt so that the Nyquist spacing is 1. We see that once the 
adjacent point spacing has reached twice the Nyquist spacing, 
the entropy has transitioned from the logarithmic scaling law 
to a volume law. The infrared scale is uj/Q — 10“^°°. 


points closer together, we find that the entropy does not 
depend very sensitively on the spacing between points. 
The logarithmic scaling of the entanglement entropy re¬ 
mains the same, and only the subleading constant is 
modified. This is in contrast with the case of under¬ 
sampling, where we saw a transition from logarithmic to 
linear growth. 

Unfortunately, with an oversampled set of lattice 
points the above numerical calculation becomes unsta- 























14 


ble for a large number of sample points. This is because 
the matrices A and S defined in section IIIIBI become 
ill-conditioned. However, it is possible to perform these 
numerical calculations in the regimes of a small amount 
of oversampling, or for only a small number of points. 

For a small amount of oversampling, that is, with lat¬ 
tice spacing between 97%-100% of the Nyquist spacing, 
we continue to find logarithmic scaling of the entangle¬ 
ment entropy with the number of points traced out: 

S{N) = Co log N + Cl (81) 

where in each case Cq S [0.333, 0.335] and Ci = 3.25. This 
result is for an infrared scale of w/H = 10“^™. 

Although we were not able to extend far into the over- 
sampling regime, if we fix a small number of points it is 
possible to calculate the entanglement entropy of these 
points for an arbitrary amount of oversampling. Fig¬ 
ure shows the entanglement entropy for 5 sample 
points as a function of the separation between the points. 
One sees from the figure that the entanglement entropy 
reaches a plateau for spacings below the Nyquist spacing, 
with the entropy depending only weakly on the spacing 
in this regime. This shows that the logarithmic scaling 
law for the entanglement entropy is also independent of 
the spacing between the points, if the spacing is below 
the Nyquist spacing. 

In figure [TOj we also see that the entropy continues to 
increase as the spacing is increased above the Nyquist 
spacing. This is because the entropy depends on local 
correlations, and as the degrees of freedom become more 
separated, more of their correlations are with the com¬ 
plementary degrees of freedom. In the 1-f 1 field theory 
we consider, these correlations decay very slowly, so that 
the entropy continues to increase all the way up to the 
infrared cutoff scale. The entropy is also slightly peaked 
at integer multiples of the Nyquist spacing: this is re¬ 
lated to the fact that the commutators vanish between 
degrees of freedom separated by multiples of the Nyquist 
spacing. 

We can also show that the plateau in the entropy at 
small spacings is not only a feature of the ground state, 
but also occurs in the thermal state considered in sec¬ 
tion jlV^ Figure [m shows the entanglement entropy of 
five points as a function of their spacing for various tem¬ 
peratures. We see that at any temperature, when the 
points are closer than a Nyquist spacing, the resulting en¬ 
tanglement entropy depends only weakly on the spacing. 
This plateau is therefore not a fundamentally quantum- 
mechanical effect: if we perform the same calculation 
with five points for the classical thermal state, we find 
the same plateau. 

The restriction to a small number of sample points is a 
consequence of numerical instability that occurs for small 
spacing and large numbers of points. In the following 
subsection we develop an analytic method that will allow 
us to consider a larger number of points, in the limit of 
small spacing. 



Ax 


FIG. 10. Entanglement entropy of five points as a function 
of the spacing between them. The horizontal axis is scaled 
by H/tt so that the Nyquist spacing is 1. For spacings be¬ 
low the Nyqnist spacing, we see a plateaning effect indicating 
that the entanglement entropy is not sensitive to the spacing 
between the points for spacings below the Nyquist spacing. 
Above the Nyquist spacing the entropy tends to increase as 
the spacing increases. In this plot the infrared to ultraviolet 
ratio is uj/fl = 10“®. 


D. Oversampling using derivatives 


The above results suggest that even if the N sample 
points are arbitrarily close together, the corresponding 
entanglement entropy will scale as ^\og{N). However, 
these results are restricted either to a small number of 
lattice points, or to very mild oversampling. To show this 
scaling behaviour we will develop a method that allows us 
to consider a larger number of sample points in the limit 
of small spacing. In the limit where the N points are 
taken coincident, sampling at N points is equivalent to 
sampling the first N spatial derivatives at a single point. 
This is related to the result that in classical sampling 
theory that instead of sampling the field at the Nyquist 
rate, one can instead sample the field and its derivatives 
at a fraction of the Nyquist rate [33]. For example, if at 
a sample point one measures both the amplitude of the 
field and its first derivative, then the samples only need 
to be taken at half the Nyquist rate. 

First, consider a situation where we sample at just two 
nearby points. We can perform the following symplectic 
transformation on the phase space variables: 


/ Hx) 
4>{x -I- Ai) 
tt{x) 

V 7r(a: -|- Ax) 


/ -^{^{x +/Ax) + (l)[x)) \ 
Ax)-(/)(x)) 
;^(7r(x -f Ax) -I- 7r(x)) 

V + -7r(x)) y 


(82) 


Thus, we see that sampling two points in the limit where 
their separation vanishes. Ax -A 0, can equivalently be 
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viewed as sampling the field and its derivative at a single 
point (as well as the conjugate momentum and its deriva¬ 
tive at that point). This motivates performing the cal¬ 
culation of the entanglement entropy after sampling the 
field and its higher derivatives {(/), </>', 4>",..., at 

a single point, which is equivalent to sampling the field 
at N points in the limit of small spacing. This can be 
thought of as an extreme case of oversampling the field. 


FIG. 11. Entanglement entropy of five points as a function 
of the spacing between them at various temperatures. The 
horizontal axis is scaled by SI/tt so that the Nyquist spacing 
is 1. We see here a plateau in the entropy at small spac- 
ings, similar to the plateau in the ground state entropy. The 
analogous classical calculation is also shown as black dots for 
temperature T /fl = 1. At temperatures T/Q. > 1, the entropy 
behaves the same as for T/fl = 1 but shifted vertically by a 
constant ~ log(r/f2). In this plot the infrared to ultraviolet 
ratio is uj/Q, = 10“®. 


The required matrix elements for the calculation of the 
entropy are: 


0 it (n -I- to) = 1 mod 2 


(83) 


(_l)"+2"’"^log(g) ifn-bTO = 0 

* ^ Tbibn (^"^™ ~ 0 ;"+™) if n -I- TO 0, (n -I- to) = 0 mod 2 , 
0 if (n -|- to) = 1 mod 2 


{{d^rm ■ {d,rm) = { (-i) 


(84) 


and 


{{dxTT^ix) ■ {dx)""TT{x)) = 


(- 1 )^ 


J_ 1 

277 n+m+2 


^i^n+Tn+2 _ ^n+m+2'^ H {n + m) = 0 mod 2 
if (n -|- to) = 1 mod 2 


(85) 


Figure]^ shows the entropy as a function of the num¬ 
ber of points in the limit of small spacing. For reference 
we also show the entropy for Nyquist spacing. Both re¬ 
sults are in agreement with the curve l/31og(fV) + c, but 
with different constants. 

This is consistent with the interpretation that both 
Nyquist spacing and the small spacing limit calculate the 
entropy of an interval of the same length, but regulated 
in a slightly different way. This different regularization in 
the small spacing limit slightly disentangles the degrees 
of freedom very close to the entangling surface. This re¬ 
duces the constant coefficient in the entropy, while keep¬ 


ing the leading logarithmic behaviour. We will give some 
evidence for this interpretation in section |Vj 


V. LOCALIZATION OF DEGREES OF 
FREEDOM 

We have shown how the entanglement entropy of a set 
of samples of a bandlimited quantum field depends on the 
number of samples taken and on the spacing between 
samples. For spacings larger than the Nyquist spacing 
the entanglement grows linearly with the number of sam- 
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FIG. 12. Entanglement entropy dependence of number of 
derivatives traced out at a single point. Fitted curve to 
Nyquist spaced points is S{N) = 0.334 log(A^) + 3.25, where 
N is the number of sampling points. Fitted cnrve to deriva¬ 
tive sampling points is S{N) — | \og{N) + 3.14, where N — 1 
is the number of derivatives sampled (which corresponds to N 
sampled points). We see that both cnrves differ by a constant 
« 0.11. The infrared to ultraviolet ratio is uj/Q, = 10“^°°. 


pies. At the Nyquist spacing there is a sharp transition. 
For samples spaced at a Nyquist spacing or closer, the en¬ 
tanglement entropy grows logarithmically with the num¬ 
ber of samples, with a coefficient that is independent of 
the spacing. We will now show how this is a consequence 
of each mode of the field occupying an incompressible 
volume of space, which we will call a Planck volume. 

Since the degrees of freedom on all sample points of a 
Nyquist lattice is equivalent to the entire field, one can 
interpret each degree of freedom on the Nyquist lattice as 
representing information roughly contained in a Planck 
volume centered at the sample point. Thus, we can in¬ 
terpret sampling a set of points on a Nyquist lattice as 
sampling the corresponding interval of the field. If we 
sample finitely many contiguous degrees of freedom on 
a Nyquist lattice, the resulting entanglement entropy is 
generated by the correlations cut across the boundary of 
the interval. Thus, the entanglement entropy varies loga¬ 
rithmically with the number of points sampled, N, since 
they describe an interval of length L = Ntt/^1. 

If the sample points we take are spaced farther than a 
Nyquist spacing, the corresponding Planck volumes cen¬ 
tered at these sample points will be disjoint. Therefore, 
the local correlations around each point which contribute 
to the entanglement entropy will mostly be independent. 
Therefore, the entanglement entropy for points which are 
well-separated will scale linearly with the number of sam¬ 
pled degrees of freedom, as we demonstrated above. 

The most interesting case is the case of oversampling. 
When the sampled degrees of freedom are pushed very 
close together, we find that the entropy varies logarith¬ 


mically with the number of degrees of freedom sampled. 
As the degrees of freedom are taken together one might 
expect to probe only a single Planck volume, and hence 
find a smaller entropy. If the N points are almost on 
top of one another, the correlations at the boundary of 
the overlapping Planck volumes would generate entan¬ 
glement entropy roughly equivalent to the entanglement 
entropy of a single sampled degree of freedom, rather 
than N degrees of freedom. However, from the results of 
the entropy calculation, it seems that the N degrees of 
freedom correspond to a larger, effective volume the size 
of the original volume described by the N contiguous 
Nyquist-spaced points before they were pushed together. 

We can estimate the size of this effective volume as 
follows. First, we take a set of N sample points with 
some equidistant spacing, which we shall denote subsys¬ 
tem Atv. Now we include another sample point which 
we initially place very far from these N sample points, 
which we shall denote subsystem P. This point will be 
used as a probe to localize the points in subsystem A^r, 
using the mutual information 

I{AN:P) = S{AN) + SiP)-S{AN,P), ( 86 ) 

where S{A]\f), S{P) are the entanglement entropies of the 
subsystems An and P (respectively), and S{An,P) is 
the entanglement entropy of the combined system. While 
the probe point P remains far away, the mutual informa¬ 
tion will be small because the correlations between An 
and P are small. When P approaches the system of N 
points, the mutual information will increase as the two 
subsystems become more correlated. We can use the in¬ 
crease of the mutual information I {An '■ P) as an indi¬ 
cator of the extent of the subsystem An- In particular, 
we can map the boundary of An for various values of the 
spacing between the N points. If the points spaced far 
below the Nyquist spacing indeed describe an effective 
volume of ~ Planck volumes in size, then we should 
see the boundary of the subsystem at the edge of this 
region. 

In figure we see this is the case for a subsystem 
of A = 5 points. When the points in the interval are 
spaced farther than a Nyquist spacing (in the upper re¬ 
gion of the graph), we can clearly identify the region of 
space that they occupy, which is centered at each sam¬ 
ple point and on the order of a Planck volume in size. 
When the spacing approaches the Nyquist spacing, the 
points begin to occupy a single interval of roughly A = 5 
Planck volumes in size. As the spacing decreases below 
the Nyquist spacing, the size of this interval ceases to 
decrease, indicating that the sample points describe the 
same volume of space regardless of their positions below 
the Nyquist spacing. Therefore, the sample points each 
describe an independent volume of space of Planckian 
size. In this sense, the Planck volumes described by the 
individual degrees of freedom are incompressible. 

This behavior shows how the field is able to have a 
finite information density, while still allowing arbitrarily 
closely spaced probes. If we attempt to sample the field 
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FIG. 13. Mutual information between subsystem Am and probe point P. The spatial axes are scaled so that the Nyquist 
spacing is 1. The variable x denotes the position of the probe point P, with x = 0 at the centre of the subsystem Am- For 
spacings much larger than the Nyquist spacing, the degrees of freedom occupy independent intervals of the order of the Nyquist 
spacing in size. For spacings at or below the Nyquist spacing, the points in the interval describe a fixed volume of size N. 


at two points closer than a Nyquist spacing, we are really 
only probing the degrees of freedom in a larger region 
of space centered around these samples, whose volume 
is determined by the number of samples rather than by 
their spacing. 


VI. INFRARED BEHAVIOUR 

Above we have examined the effect of the ultraviolet 
cutoff on the calculation of the entanglement entropy of a 
bandlimited quantum field. Here we also briefly examine 
the infrared behaviour of the entanglement entropy. Refs. 
[43l |44] calculated the entanglement entropy of a free 
scalar field in 1+1 dimensions with an infrared cutoff. 
They find that the entanglement entropy of a subset of 
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N oscillators takes the form 

(i))-log,, (87) 

where cq is a constant and w is the infrared cutoff. Here 
the entropy diverges as a double logarithm with the in¬ 
frared cutoff. This infrared divergence is a special fea¬ 
ture of H-l dimensions and arises due to a build-up of 
long range correlations in the field. We saw in section 
|TT] that the (j)-(j) correlations decay very slowly with the 
point separation, and are only suppressed on the order of 
the longest wavelength 27r/a;. If we remove the infrared 
cutoff, we allow for infinitely long wavelengths, and thus 
divergent correlations. 

In our above numerical results for calculating the en¬ 
tanglement entropy of a quantum field on a Nyquist lat¬ 
tice, we perform a numerical fit of the data to determine 
the coefficients of the equation 


arises because the diagonal elements of the (j)-(j} correla¬ 
tion matrix diverge logarithmically in the infrared as 



One may wonder how the infrared behaviour of the en¬ 
tanglement entropy changes with a different implementa¬ 
tion of the infrared cutoff. In particular, we will examine 
how the entropy behaves when the field modes are dis¬ 
crete. Consider a scalar field on an interval [0, L] with 
periodic boundary conditions (up to a phase), 

,^(0) = e*“(/>(L), (96) 

where a G [0,27r). Expanding the field in a spatial 
Fourier series. 


S' = Co log(A^) + 2 log ( log 


S(fV) = Cl log iV -I- C 2 . (88) 

This numerical fit was also performed for various values 
of the infrared cutoff, uj/fl € [I0“®°, 10“^™]. We find 
that the leading order behaviour of C 2 in this range is 

02 = (log (log (0)). (89) 

The coefficient C 2 agrees with the result of Ref. [H] (up 
to a constant). Furthermore, in the limit u; ^ H, the 
infrared cutoff oj has the same effect on the entropy as a 
small mass [15] . 

It is possible to anticipate this result analytically by 
considering tracing out a system of one oscillator. From 
the expressions found above for the two-point functions, 
we find the diagonal terms are 

l(/)(x),7r(x)] = i- —(90) 

TT 

(02(a:)) = ^ log , and (91) 

(7r^(x)} = ^ (fl^-uj^) . (92) 

The symplectic eigenvalues of the covariance matrix are 

Thus, the entanglement entropy is 
S' = (A+ -k ^) log(A+ + ^) - (A+ - ^) log(A+ - i), (94) 


Hx) = 

k 


Akx 


y/uJkL 


(f’k 


(97) 


where uik ■= |fc|, one obtains a discrete set of momentum 
modes {kn = }ngz- The ultraviolet cutoff will be 

imposed by restricting the set of allowable n to {n G 
Z : ujk^ = \kn\ < H}. We see that for a > 0 there 
will be no zero mode, so the (fi-cj) correlations are finite, 
but there will be an infrared divergence as a —>■ 0. It 
is straightforward to find that the Hamiltonian for this 
Klein-Gordon field can be expressed as 

H ■ (98) 


Expanding in terms of the usual creation and annihilation 
operators. 


ak 


1 

71 


{4>k + jTTfc) , 


one finds 


(99) 


[0((r),7r(cr')] = ^^e*""("-"7 (100) 

n 

AiAAx')) = and (101) 

2L bjh 

n ^ 

{Ax)Ax')) = ^ ( 102 ) 


which to leading order in the limit w/H —>■ 0 is 

S~ll„g(l„g(^)). (95) 

We see directly that the outer logarithm is due to the 
entanglement entropy formula, and the inner logarithm 


These matrix elements are easy to calculate numerically 
since the ultraviolet cutoff imposes a finite number of n 
to sum over. 

Similar to the calculation performed above, we shall 
calculate the entanglement entropy associated with trac¬ 
ing out the field oscillator at a single point. The relevant 
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diagonal elements of the correlation matrices are: 
2K 


(a:),7r(a;)] = i- 


L 


(103) 


1 


1 r 




a 


K- ^ + 1 
27T 

a 
27T 




’M) = 


ttK^ 


(104) 

(105) 


where K := nL/27r, and ipix) is the digamma function. 
We see that now instead of a logarithmic infrared diver¬ 
gence of {(j)‘^{x)), with a discrete set of modes the correla¬ 
tion function diverges as 1/a as a —>■ 0. In this limit, the 
leading order behaviour of the symplectic eigenvalues is 


(106) 



The corresponding entanglement entropy is 




(107) 


Numerical calculations were performed with multiple 
points traced out, and the observed behaviour in the limit 
a —)■ 0 agrees with that in equation (107). 


For our ultraviolet and infrared bandlimited quantum 
field on the real line, the maximum allowable wavelength 
is Xmax = 2 tt(uj. Foi' periodic boundary conditions, the 
maximum wavelength is Xmax = T/27ra. Thus, we can 
compare and summarise the infrared behaviour of the 
entanglement entropy as 


• For a continuous set of modes: 

S' ~ ^ log (log (flA 

max )) as A 

n 

• For a discrete set of modes: 

S-^log^^l asA„ 


00 . 


(108) 


(109) 


In the case of boundary conditions that are periodic up 
to a phase, the entropy diverges more quickly. This is 
plausible since this case is a description of a field on a 
circle, and any two points on this circle coupled to one 
another an arbitrary number of times in the limit Xmax 
00 . 


VII. OUTLOOK 

Shannon sampling, which in information theory estab¬ 
lishes the equivalence of continuous and discrete informa¬ 
tion, also yields a simple model for how a natural ultra¬ 
violet cutoff at the Planck scale might manifest itself in 


QFT. Namely, when approaching the Planck scale from 
lower energies, where QFT is still valid, quantum fields 
are modeled as being bandlimited. 

Within this model, we probed the localization of the 
quantum field’s degrees of freedom by tracking the be¬ 
haviour of the vacuum entanglement entropy. We found 
that the degrees of freedom of the quantum field, i.e., 
the local field oscillators, are non-local but can be local¬ 
ized down to the Planck-scale. In fact, we found that the 
local degrees of freedom occupy incompressible Planck- 
scale volumes, in the sense that N degrees of freedom 
always describe N Planck volumes regardless of the dis¬ 
tance between them. 

The tools that we outlined in this paper are applicable 
also to the study of bandlimited quantum fields in higher 
dimensions, and to fermionic fields. 

For the case of higher dimensions, we conjecture that 
we will continue to see incompressibility of the field de¬ 
grees of freedom, i.e., a plateau in the entanglement en¬ 
tropy for spacings below the Nyquist spacing. A system 
of N lattice points should continue to probe a region of 
N Planck volumes even if they are sub-Nyquist spaced. 
In higher dimensions it will be interesting to see how the 
shape of this region changes with the relative location 
of the sample points. Moreover, at larger than Nyquist 
spacings, we expect to find an infrared plateau. This is 
because the vacuum correlations decay over shorter dis¬ 
tances in higher dimensions and do, therefore, not exhibit 
the infrared divergences that soften the area scaling to a 
logarithmic scaling in 1-1-1 dimensions. 

Our approach to studying localization and vacuum en¬ 
tanglement should also work for quantum fields on curved 
spacetimes. To this end, the generalization of Shannon 
sampling theory to curved spaces developed in [3S] can 
be used. It will be interesting, in particular, to determine 
how the localizability of the fundamental field oscillator 
degrees of freedom is affected by curvature and horizons, 
and in particular, how this may affect Hawking radiation. 
A similar hard momentum cutoff was considered in [45) . 
and it should be interesting to apply our new methods 
to that case. 

A further direction of great interest would be to study 
the entanglement entropy and the localization of degrees 
of freedom of fields with a fully covariant ultraviolet cut¬ 
off. In this case, the bandlimitation would be imposed 
by cutting off the spectrum of a covariant Laplacian 
(on Riemannian manifolds) or a d’Alembert operator (on 
Lorentzian manifolds) |46j . To this end, it may be use¬ 
ful to express the entropy purely in terms of spacetime 
correlation functions [42], i.e., without reference to the 
canonical formalism. We remark that it has been argued 
that if the entanglement entropy is finite and obeys a 
form of the Clausius relation, then Einstein’s equation 
emerges as a thermodynamic equation of state [13 |48] . 
Thus having a regulator that is both covariant and cuts 
off the entanglement entropy could be a key step toward 
quantum gravity. 

It should also be very interesting to study the inter- 
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action of Unruh-DeWitt detectors with the bandlimited 
quantum fields. Unruh deWitt detectors would describe 
explicit means to sample the quantum fields, and this 
could, therefore, constitute a further step towards the 
quantization of sampling theory. An interesting question 
that then arises is how the spatial profile of an Unruh- 
DeWitt detector interacts with the spatial profile of the 
degrees of freedom. Another key question that arises 
with the quantization of sampling theory follows from the 
fact that, in classical sampling theory, the reconstruction 
of a function from a lattice with non-equidistant sample 
points is more sensitive to noise in the sample measure¬ 
ments than reconstruction from a lattice with equidis¬ 
tant sample points. This, therefore, raises the interesting 
question if or to what extent quantum noise may inter¬ 
fere with the ability to use significantly irregularly-spaced 
lattices. 

One of our central results is that the entanglement en¬ 
tropy does not decrease significantly as the spacing is 
decreased to below the Nyquist spacing. In classical sam¬ 
pling theory too a curious phenomenon is known to arise 
as the spacing is decreased to below the Nyquist spacing, 
namely the phenomenon of superoscillations. 

Superoscillations in a bandlimited function are a finite 
set of oscillations that oscillate faster than the highest 
Fourier component in the signal [SU . Superoscil¬ 

lations are difficult to generate and rarely occur naturally 
but they have become an active held of investigation, also 
in engineering, because of their potential for superreso¬ 
lution. An open problem in this held is to quantify the 
prevalence of superoscillations in random signals. In our 
analogous classical calculation, we obtained probability 
distributions for observing particular held values on an 
arbitrary sampling lattice. These distributions could be 
used to calculate the probability for hnding the held in 
any particular superoscillatory conhguration. It might be 
possible to hnd in this way a practical measure of classical 
or quantum held huctuations that is sensitive specihcally 
to the occurrence of superoscillatory behaviour. 

A key question that we addressed in this paper has 
been the question of how a bounded continuous region 
can support only a hnite number of degrees of freedom. 
In our model for a natural ultraviolet cutoff, this is ac¬ 
complished by the quantum held degrees of freedom at 
distinct spacetime points not being independent. As a 
consequence, probing with too many operators in a small 
region simply causes one to probe a larger region of space. 
This means that the degrees of freedom of the theory are 
encoded redundantly on a single spatial slice. 

In fact it has recently been argued that, similarly, sig- 
nihcant redundancy is present in the bulk held theory in 
AdS/CFT [53]. The redundancy implied by holography 
appears to be more extreme, however, since the number 
of degrees of freedom is bounded by area rather than vol¬ 
ume. To obtain this more drastic reduction, the nontriv¬ 
ial commutation relations coming from gravitation inter¬ 
action may provide some clues [25] . One may also take 
a cue from holography by trying to redundantly encode 


the degrees of freedom of a quantum held in a space of 
higher dimension. 

Finally, we note that the hyperbolic geometry of anti- 
de Sitter space suggests that an encoding based on 
wavelets may have some connection to holography IMl- 
ETj. Wavelet theory and Shannon sampling theory are 
closely related and the quantization of wavelet theory 
should be possible along the lines that we developed here. 
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Appendix A: Functional analytic structure of 
sampling theory 

Here we will briefly review the functional analytic 
structure of sampling theory, as first shown in [19] . The 
functional analytic view is powerful because it ultimately 
reduces sampling theory to the simple fact that when 
a Hilbert space vector is known in one basis then it is 
known in all bases. A key point here is the distinction 
between symmetric and self-adjoint operators. By the 
spectral theorem, each self-adjoint operator possesses a 
unique diagonalization, whereas a simple symmetric op¬ 
erator does not, see, e.g., [5S]. A certain type of sim¬ 
ple symmetric operator, however, possesses a family of 
self-adjoint extensions, and their diagonalizations and 
spectra provide the sampling lattices of sampling the¬ 
ory. (The functional analytic view therefore also allows 
one to work with varying Nyquist rates, which can occur 
in signal processing as well as in curved spacetimes.) 

Indeed, when acting upon the space of functions ban¬ 
dlimited by some maximum frequency D, the usual po¬ 
sition operator from first quantization, A, is symmetric 
but not self-adjoint. To see this, recall that for an oper¬ 
ator X to be symmetric, it must obey 

= (<^|AV^) (Al) 

for all \(j)),\ilj) € D{X), where D{X) is the domain of 
X. Equivalently, an operator is symmetric if and only 
if its expectation value is real for all vectors in its do¬ 
main: {(j)\X(j)) G TZ Vcj) G D{X). For an operator 
to be self-adjoint, it must be symmetric and its do¬ 
main must coincide with the domain of its adjoint, i.e., 
{X(l)\tP) = {(l)\X'ilj) V|()>), |^/>) G D{X) and D{X) = D{X^). 

Indeed, the position operator X with domain restricted 
to the space of bandlimited functions is symmetric, but 
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the domain of its adjoint is a larger space of functions. 
This is most easily seen in the momentum eigenbasis. 
In this basis, the space of physical wavefunctions, and 
therefore the domain of X, is the space of functions 

D{x) = {(^ e L‘^[-n,n]\(i) e 

(I)' G = 0 } (A2) 

where AC\—Vl^ fl] denotes the space of absolutely contin¬ 
uous functions. Notice that D{X) contains only func¬ 
tions which obey Dirichlet boundary conditions. In con¬ 
trast, the domain of X^ is the larger function space ob¬ 
tained by not imposing any boundary conditions. This 
can be seen from the definition of : 

= j dk ^*{k) 

=i (j)* 

+ j dk {k)^{k). 

We see from the last line that as long as |'0) is in the do¬ 
main of X (thus obeys Dirichlet boundary conditions on 
the interval [—D, D] in Fourier space), then we can define 
the adjoint of X^ as i{d/dk) in the momentum represen¬ 
tation. However, since we did not need to impose bound¬ 
ary conditions on \(j)), we have that 11(A) C D{X^). 

Crucially now, it is possible to extend the domain of 
X so that the extension is a self-adjoint operator. This 
operator is called a self-adjoint extension of X (e.g., see 
[SSllSS]). In our situation, we can perform this extension 
by enlarging the domain of X to include functions with 
periodic boundary conditions up to a phase on the in¬ 
terval [—D,D] in Fourier space, i.e., = e“*“'0(—fl). 

We shall denote the self-adjoint extension corresponding 
to the particular phase as thus the family 

of self-adjoint extensions is parametrized by a S [0,27r). 
Now we will show that the operators are indeed 

self-adjoint. The definition of the adjoint of gives 


Now, since the operators X^°‘'> are self-adjoint, they 
have spectral decompositions. The spectrum of the oper¬ 
ator X^°‘'> (for fixed a) is discrete, spec(A(“)) = := 

and describes a one-dimensional lattice. The 

corresponding eigenvectors, {|a;l“^)}„gz, are represented 
in the momentum eigenbasis as 

(fc|x(“)) = (A5) 

These eigenvectors are orthogonal and admit a resolution 
of the identity: 


^|a:(“))(4“)|=l (A6) 


. Crucially, the position eigenvectors from different self- 
adjoint extensions are not orthogonal: 


( 4 . 


\444) =sinc [(4“^ 


(A7) 


Note that the union of the spectra of the entire family of 
self-adjoint extensions provides a covering of IR. There¬ 
fore, it is possible to construct an overcomplete contin¬ 
uum basis by taking the union of eigenbases of the family 
of self-adjoint extensions, i.e., \x) := \x4'^) x = 

„(“)_ 277 - 71 - 1 - 0 : 
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It is then simple to write down the Shannon sampling 
theorem for a bandlimited function "0: 


(xlV') 

(AS) 

^(x|xl“))(x(“)|^) 

(A9) 



^sinc {x-x4^)n 'ipix^'’)- 

(AIO) 




Therefore, we see that the function ip is determined at 
any point x G K, from its values on one of the lattices 

{x44n- 

Also, we obtain an overcomplete resolution of identity. 




+ J dk (fc)V'(fc) 

=i (^*(D) - e“(^*(-D)) i/)(D) 

+ J dk (fc)V'(fc)- 


(A4) 


We see that the adjoint of is defined as i{d/dk) in 
the momentum basis provided that \(j)) also obeys the 
boundary condition d)(D) = Therefore we 

see D{X^°‘'>) = D{X^°'>^), and so is self-adjoint. 


n 

TT 


dx |x)(x| = 1 




(All) 


where tt/D is the density of degrees of freedom. We can 
use the resolution of identity for the continuum basis to 
show that the space of bandlimited functions has a re¬ 
producing kernel, 

(x|^) = —y dx' (x\x'){x'\'ip) (A12) 

'ipix) = J dx K{x,x')ijj{x'), (A13) 

where K{x,x') := (D/7r)(x|x') = (D/tt) sinc[(x — x')D]. 
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